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EXECUTIVE  SUMMARY 

The  question  of  how  much  error  is  present  in  a  simulated  infrared 
signature  if  reflections  are  neglected  was  studied  using  a  sensitivity 
analysis  approach.  The  purpose  of  this  study  was  to  determine  whether  it 
is  necessary  to  model  reflections  in  the  current  FSTC  method  of  simulating 
IR  signatures.  To  gain  more  insight  the  following  steps  were  performed: 

1)  A  sensitivity  analysis  approach  was  chosen  for  the  study. 

2)  An  equation  isolating  the  major  components  of  an  IR  signature  was 
derived, 

3)  The  equation  was  used  to  determine  error  present  in  a  simulation 
when  emissivity  was  estimated  to  be  always  0.95.  This  was  accom¬ 
plished  by  using  a  simulation  of  a  complicated  tank  applied  with 
actual,  measured  emissivity  that  had  specular  characteristics^ 

4)  A  weighted  average  error  figure  was  derived  and  applied  to  the 
results  of  histogramming  emissivities  for  the  simulated  tank^ 

5)  It  was  determined  that  in  some  cases,  error  could  be  as  high  as 
35%  on  an  absolute  scale.  However,  most  of  the  errors  calculated 
in  section  3.0  were  found  to  be  less  than  10%,  which  corresponds 
to  about  6*K  apparent  temperature  error  for  ambient  targets. 

Average  error  for  a  typical  target  was  found  to  be  on  the  order  of 
6*K  apparent  temperature  when  its  IR  signature  was  simulated  using  FSTC's 
current  method.  FSTC  personnel  must  judge  whether  6*K  apparent  tempera¬ 
ture  error  is  significant  before  further  action  is  planned. 
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1.0 

WHY  A  SENSITIVITY  ANALYSIS  OF  IR  SIMULATION? 

This  report  is  the  second  in  a  series  addressing  the  issue  of  how 
much  error  is  introduced  into  IR  target  simulations  by  neglecting  the 
effects  of  reflections  from  a  target's  environment.  Technical  Memo  #1  of 
this  contract  explored,  through  a  case  study  analysis,  reflections  in  the 
IR  and  the  error  caused  by  neglecting  them  in  a  simulation.  This  report 
explores  the  same  issue  through  a  sensitivity  analysis  approach. 

A  rigorous  simulation  of  IR  signatures  should  include  not  only  target 
self-emission,  but  reflections  from  its  environment  and  itself.  However, 
neglecting  reflections  reduces  the  cost  of  software  development  and  data 
preparation.  The  question  of  whether  reflections  need  to  be  included  in 
an  IR  simulation  is,  thus,  an  important  one  since  it  is  recognized  that 
target  emission  is  often  dominant  in  the  IR  wavelengths.  Whether  or  not 
the  additional  expense  of  including  reflections  is  justified  is  ultimately 
determined  by  the  required  simulation  accuracy.  An  ideal  situation  would 
be  a  simulation  that  is  "accurate  enough"  for  a  given  purpose  and  models 
only  target  emission,  in  order  to  avoid  more  expensive  development  and 
preparation.  It  is  the  intent  of  this  report  to  illustrate  the  magnitude 
of  error  that  results  from  ignoring  reflections.  The  significance  of 
these  errors  must  be  judged  by  FSTC  personnel. 

Approaches  to  Modeling  Emissivity  and  Reflections 

In  general,  surfaces  can  be  modeled  as  either  perfectly  specular, 
perfectly  diffuse,  or  somewhat  intermediate  to  these  extremes.  Diffuse 
surfaces  with  high  emissivity  (>0.95),  can  be  modeled  with  "-S’K  apparent 
temperature^  error  [3].  Depending  on  the  use  of  the  simulation,  this  may 
be  adequate  or  inadequate.  FSTC  currently  assumes  that  this  error  is 


1.  Apparent  temperature  is  the  temperature  of  a  perfect  emitter  whose 
radiance  equals  that  received  from  the  target. 
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acceptable  and  assumes  that  most  modeled  materials  are  diffuse  and  have  an 
emissivity  close  to  0.95. 

However,  for  surfaces  that  are  significantly  specular,  not  diffuse  as 
modeled  by  FSTC,  certain  viewer/target  geometries  will  result  in  a  large 
amount  of  reflections.  Neglecting  reflections  when  simulating  an  image  of 
a  specular  target  at  certain  viewer/target  geometries  may  result  in  large 
and  even  unacceptable  error.  The  question  to  be  answered  is:  How  often 
and  how  large  will  the  error  be  from  neglecting  reflections  when  simu¬ 
lating  IR  imagery? 

A  Sensitivity  Analysis  Applied  to  IR  Image  Simulation 

In  order  to  achieve  a  better  understanding  of  the  importance  of 
various  factors  in  an  IR  image,  a  sensitivity  analysis  was  used.  The 
analysis  was  performed  using  the  factors  known  to  be  important  in  an  IR 
image.  These  factors  are: 

•  emissivity 

•  surface  temperature 

•  reflection  component. 

Figure  1  depicts  the  three  components  of  an  IR  image. 

Emissivity  is  a  function  of  surface  material  properties,  wavelength 
and  viewer/target  geometry  (if  it  is  a  material  that  has  significant 
specular  behavior).  Most  military  type  paints  exhibit  some  specularity, 
i.e.,  emissivity  decreases  with  decreasing  grazing  angle  as  shown  in 
Figure  2.  This  fact  in  large  part  has  prompted  the  present  study. 

Surface  temperature  is  used  to  compute  the  emitted  radiance  leaving 
the  target.  The  relationship  between  temperature,  wavelength  and  emitted 
radiance  is  described  by  Planck's  law.  FSTC's  simulation  uses  Planck's 
law  to  compute  emitted  radiance. 

Reflection  component  is  often  characterized  by  an  apparent  source 
temperature  and  the  true  emissivity  of  the  target.  Thus,  a  clear  night 
sky  is  often  said  to  be  at  an  apparent  temperature  of,  say,  240*K,  and  its 
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FIGURE  1.  THE  THREE  COMPONENTS  OF  AN  IR  SIMULATION 
STUDIED  IN  THIS  REPORT. 
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FIGURE  2.  FACTORS  THAT  AFFECT  THE  VALUE  OF  EMISSIVITY. 
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radiance  is  obtained  from  Planck’s  law,  similarly  to  surface  temperature. 
This  radiance  is  then  reflected  off  the  target  in  an  amount  determined  by 
the  latter's  emi ssi vi ty/refl ecti vi ty  nroperties.  For  the  predominantly 
specular  target  considered  in  this  report,  only  a  single  source  direction 
need  be  considered.  In  the  general  case,  the  reflection  component  of  a 
target  signature  includes  contributions  from  an  entire  hemisphere  of 
source  temperatures  (cf.,  equation  A-7  in  the  Appendix). 

The  purpose  of  the  sensitivity  analysis  is  to  determine  the  magni¬ 
tudes  of  error  that  could  be  expected  in  the  current  simulation  method 
used  by  FSTC.  The  current  method  assumes: 

1)  Uniform  emissivity  (i.e.,  independent  of  view  angle),  typically 
0.95  for  8-14  fim  simulations  and  0.6  for  3-5  /xm  simulations. 

2)  Reflections  of  any  sort  are  an  insignificant  component  of  the 
total  signature. 

Both  of  these  assumptions  are  open  to  question.  The  current  study  is 
intended  to  illuminate  exactly  how  erroneous,  or  how  accurate,  these 
assumptions  can  be. 
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2.0 

COMPUTING  SENSITIVITY  TO  ERROR 

An  equation  that  relates  the  error  in  final  radiance  given  that  etnis- 
sivity,  surface  temperature,  and  the  reflection  term  have  been  incorrectly 
estimated  is: 


AL 


(e  -  e)  L(T^)  +  e 


31 

< 

3T 

^T, 

AT  -  (1  -  e)  L(Tg) 


(1) 


where: 

e  =  true  emissivity  of  the  target 
=  f(AX.fl), 

AX  =  integrated  wavelength  band, 

Q  -  grazing  angle  between  viewer  and  target, 

A 

e  =  estimated  emissivity  of  the  target, 

L()  =  Planck  operator  that  computes  radiance  based  on  temperature  and 
radiance, 

=  apparent  temperature  of  source  of  reflection, 

Tj  =  actual  temperature  of  target, 

AT  =  (T  -  T)  or  difference  between  estimated  temperature  and  actual 
value. 

The  derivation  of  equation  1  is  presented  in  Appendix  A.  Equation  1 
is  used  as  the  basis  of  the  following  sensitivity  analysis. 

Each  of  the  th»'ee  terms  in  equation  1  can  be  interpreted  as  being  the 
error  contribution  for  each  of  the  three  factors  in  an  IR  image  identified 
in  section  1.0.  The  error  due  to  incorrectly  assigning  emissivity  is: 

(e  -  e)  L(T^) 

The  error  due  to  incorrectly  assigning  surface  temperature  is: 
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and  the  error  due  to  reflections  is: 

(1-e)  L(T^) 

2.1  REDUCING  THE  RADIANCE  EQUATION  TO  AN  ERROR  EQUATION 

Absolute  radiance  error  in  terms  of  physical  units  lends  little  in¬ 
sight  into  the  question  of  how  “big"  an  error  is  actually  present.  One 
way  of  introducing  relativity  is  to  normalize  equation  I  in  order  to  pro¬ 
duce  a  unitless  quantity  on  the  left  side.  The  question  is,  normalize 
with  what  quantity?  An  obvious  answer  would  be  to  normalize  with  L(Ty)  to 
produce: 


&L 


A 


[3t1TtJ 


AT  -  (1  -  e) 


L{T  ) 

rriyT 


(2) 


where  now,  AL/i(Tj)  =1.0  would  be  interpreted  as  100%  error. 

2.2  SELECTING  FOUR  CASES  OF  REFLECTANCE  STRENGTH 

Equation  2  is  difficult  to  understand  because  of  the  complex  manner 
in  which  the  various  terms  can  interact  with  one  another.  One  way  of 
reducing  this  difficulty  is  to  select  one  or  two  of  the  variables  in  equa¬ 
tion  2  to  be  looked  at  discretely,  that  is  select  a  few  cases  and  then 
vary  the  other  variables  continuously.  As  mentioned  earlier,  Technical 
Memo  #1  pointed  to  emissivity  as  a  driving  factor  in  error.  Therefore, 
L(T^)/L{T.j.)  will  be  set  at  a  few  discrete  values  and  emissivity  and  other 
terms  will  be  varied.  This  will  be  used  to  establish  a  trend  in  error 
calculations  as  a  result  of  view  angle  dependent  emissivity. 

What  values  of  L{T^)/L(T.j.)  to  choose?  Two  limiting  values,  such  as 
L(Tg)/L(T.j.)  =  1  and  0,  are  useful  in  order  to  gain  insight  into  error. 
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Two  other  intermediate  values  will  fill  out  the  solution  space  suffi¬ 
ciently  to  understand  how  the  error  will  vary.  For  ease  of  discussion, 
1  et 


's' 

=  Rp  {Reflectance  ratio} 

When  =  1,  the  implication  is  that  the  source  of  reflection  =  T-j.,  the 
target  temperature.  When  R^  =  0,  the  implication  is  that  or 

that  the  source  of  reflection  is  much  less  than  the  target.  This  would  be 
typical  of  a  hot  exhaust  stack.  For  the  in  between  values,  it's  more 
difficult  to  relate  some  actual  scenarios  to  some  typical  values.  If  R^  = 
0.5  is  chosen,  one  situation  that  can  cause  this  is  T„  =  275  and  Tt  =  320 
(among  others).  This  could  be  typical  of  a  warm  vehicle  sitting  in  an 
open  field  during  the  winter.  If  R^  =  0.667,  a  situation  that  could  cause 
this  is  Tg  =  300  and  T^  =  320,  which  could  be  typical  of  a  warm  vehicle 
sitting  in  an  urban  setting.  Other  interpretations  could  be  applied. 

Table  1  summarizes  the  reflectance  ratios  chosen. 


TABLE  1 

Reflectance  ratios. 

L(T^)/L{Tt) ,  chosen  and  their 

applied 

meanings 

S  1 

R 

r 

=  L(T^)/L(T^) 

1  meaning 

R^ 

=  1.0 

Target  Self-illumination  or 

1 

source  at  same  temperature 

as  target 

R 

r 

=  .667 

Urban  Setting 

R 

=  0.5 

Open  field,  winter  time 

r 

setting  under  a  clear  sky 

R 

=  0.0 

Source  «  target  (i.e.. 

exhaust  stack) 

^RIM 


2.3  WHEN  TEMPERATURE  IS  INFERRED  EXACTLY 

Equation  2  presented  above  is  still  cumbersome  to  study  because  of 
the  term  containing  (9L/8T).  A  simpler  form  of  equation  2  to  work  with  is 
one  that  contains  error  terms  for  emissivity  and  reflection  or'y. 
Justification  for  dropping  the  second  term  in  equation  2  comes  from  Tech¬ 
nical  Memo  #1  of  March  1988,  where  it  was  implied  that  the  current  method 
FSTC  uses  to  infer  temperature  results  in  close  approximations  of  surface 
temperature.  The  greatest  error  introduced  was  found  to  be  the  error  in 
emissivity  when  a  signature  was  being  calculated  from  a  different  view 
angle  and  wavelength  band  than  the  one  at  which  the  original  data  were 
collected.  Further  analysis  will  be  presented  here. 

From  Appendix  A,  an  equation  describing  AT,  the  error  in  inferred 
temperature,  can  be  written  as 


AT  = 


e  L(T^)  +  (1  -  €)L(T5)} 


(3) 


where: 

T^  =  true  target  temperatures, 

e  =  estimated  emissivity, 

e  =  true  emissivity, 

T^  =  reflection  source  temperature, 

L~^()  =  temperature  inferred  from  a  given  radiance, 
where  all  of  the  above  quantities  pertain  to  the  waveband/view  geometry  at 
which  the  target  radiance  was  actually  measured  (cf..  Appendix  A). 

For  the  present  study,  ?  =  0.95,  and  since  0%  error  is  desired,  AT  =  0  as 
well.  Then 


Tt  =  L 


-1 


1.05  {  eL(T^)  Ml  -  e)  L(T^)}  ] 


(4) 
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Or,  rewriting  in  terms  of  radiance 

L(T^)  =  1.05  {eL(T^)  +  (1  -  e)  L(T^)}  (5) 

and  then  dividing  by  L(T-j.) 


where,  L(T^)/L(Ty)  =  R^,  the  reflectance  ratio  presented  earlier. 
Choosing  the  same  values  of  R^,  namely  R^  =  1,  0.5,  0.667  and  0,  allows 
equation  6  to  have  a  single  variable,  e.  Table  2  shows  the  emissivities 
that  result  in  0%  error  in  inferred  temperature  when  the  emissivity  is 
estimated  to  be  0.95  for  the  reflectance  ratios  chosen. 


Emissivities 

temperature 

TABLE  2 

that  result  in  0%  error  in  inferred 
vs  Reflection  ratio  R^  =  L(T^) /L(T-|.) 

R 

e 

r 

1 

*  impossible 

.667 

.85 

.5 

.9 

0 

.95 

Table  2  shows  that  when  the  reflection  source  and  target  temperature 
are  the  same,  no  real  emissivity  combined  with  the  assumed  emissivity  of 
0.95,  can  result  in  0%  error,  which  makes  intuitive  sense.  When  the 
reflection  source  temperature  is  lower  than  the  target  temperature  (as  is 
usually  the  case),  the  actual  target  emissivity  must  be  relatively  high 
(>0.85)  in  order  to  have  small  errors  in  inferred  temperature. 
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The  importance  of  this  analysis  shows  that  the  actual  target  emissiv- 
ities  must  be  relatively  high  for  there  to  be  0%  error  in  inferred  tem¬ 
perature  when  a  high  estimated  emissivity  is  used.  Naturally,  this  is  in 
direct  opposition  to  the  basic  assumption  of  this  report,  i.e.,  that  emis¬ 
sivity  at  some  view  angles  may  actually  be  low  due  to  the  specular  nature 
of  real  world  target  surfaces.  However,  one  must  keep  in  mind  that  we  are 
presently  considering  two  different  emissivities.  One  is  the  emissivity 
of  the  target  at  the  waveband/geometry  from  which  actual  target  radiance 
is  measured.  This  is  denoted  in  Appendix  A,  and  is  the  emissivity 
relevant  to  the  error  in  inferred  target  temperature  AT.  The  second  emis¬ 
sivity,  denoted  £2  ^^6  appendix,  is  at  the  waveband/geometry  of  the  IR 

simulation.  One  may  reasonably  hypothesize  a  situation  in  which  is 
high  (i.e.,  target  radiance  measured  at  near  normal  incidence)  but  £2  is 
low  (e.g.,  simulation  view  angle  near  grazing).  Indeed,  the  derivation  of 
equation  (1)  assumed  the  error  in  inferred  temperature  was  not  large;  this 
can  only  occur  if  is  reasonably  high. 

A  given  error  in  inferred  temperature  AT  represents  an  essentially 
fixed  amount  of  error  in  the  simulated  radiance.  Thus,  in  order  to  para¬ 
metrically  exercise  our  error  equation  (2),  we  may  fix  the  term  involving 
AT  to  any  desired  level.  For  the  sake  of  simplicity,  we  choose  this  level 
to  be  zero  in  the  remainder  of  this  report.  This  permits  the  impact  of 
the  remaining  terms  to  be  more  easily  understood. 

2.4  THE  REDUCED  EQUATION  TO  BE  STUDIED 

With  the  dropping  of  the  second  term  that  describes  the  error  due  to 
error  in  inferred  temperature,  the  equation  describing  error  now  becomes 


AL 

'-ri) 


=[£-£)-(!-£)  R^ 


(7) 


n 


Equation  7  will  be  used  in  the  rest  of  this  report  in  order  to  study  error 
due  to  neglecting  reflections  and  misestimating  emissivity. 

An  important  point  about  equation  3  is  that  if  the  estimated 
emissivity,  e,  is  greater  than  the  true  emissivity,  e,  and  close  to  1,  and 
the  reflectance  ratio  is  non-zero,  the  two  terms  in  equation  3  tend  to  be 
compensating.  That  is,  if 

A  ^ 

e  >  e  and  e  ~  1  and  ~  1  , 

r 


then 

[e  -  e]  -  (1  -  e)  ~  0 

This  immediately  ensures  that  error  will  be  low  as  long  as: 

1)  The  reflection  source  is  close  to  the  temperature  of  the  target 

(R^  ~  1). 

2)  The  dropped  term  concerning  inferred  temperature  is  indeed  close 
to  0. 
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3.0 

APPLYING  THE  ERROR  EQUATION  TO  A  REALISTIC  TARGET 

3.1  INTRODUCTION 

In  section  2,  a  means  of  studying  error  has  been  derived  in  equation 
form,  equation  7.  This  section  discusses  exactly  how  to  use  it  to  gain 
insight  into  sources  and  magnitudes  of  error.  This  was  partly  accom¬ 
plished  in  section  2.3  by  selecting  four  values  of  reflectance  ratio,  R^, 
in  order  to  set  limits  based  on  the  given  underlying  assumptions.  The 
task  discussed  here  is  to  further  develop  a  means  of  achieving  insight 
into  the  problem  of  reflected  term  error  iji  real istic  situations. 

3.2  A  PROCEDURE  FOR  LOOKING  AT  REALISTIC  SITUATIONS 

A  realistic  situation  would  be  to  apply  the  reduced  form  of  equation 
7  to  an  actual  target  and  derive  some  meaningful  statistics.  One  of  the 
primary  concerns  at  ERIM  is  that  since  some  target  surfaces  can  be  highly 
specular,  simulations  of  targets  at  some  view  angles  could  produce  large 
amounts  of  error  if  reflections  are  not  considered. 

Calculating  equation  7  using  measurements  from  a  foreign  target  would 
be  the  most  direct  means  of  studying  the  effects  of  neglecting  reflectance 
error.  Since  radiance  measurements  of  an  actual  tank  would  be  costly, 
resort  was  made  to  using  a  computer  simulation. 

The  GIFT  code  for  describing  and  raytracing  complicated  targets  was 
used  for  the  simulation.  A  solid  model  of  a  T72  developed  by  BRL  was 
"painted"  with  Mil-E-46096  olive  drab  paint.  The  paint's  specular  charac¬ 
teristics  had  previously  been  measured  by  ERIM  [1].  The  T72  was  then  ray- 
traced,  using  the  GIFT  code,  and  histograms  of  emissivities  were  con¬ 
structed  for  various  azimuth  and  elevation  angles.  These  histograms  were 
then  used  to  derive  error  statistics  that  lend  insight  into  how  much  error 
could  be  expected  and/or  tolerated  by  applying  the  derived  emissivities  to 
equation  7.  The  histograms  were  constructed  by  grouping  emissivities 
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computed  through  the  raytrace  into  four  categories.  These  categories  are 
shown  in  Table  3. 


Bins  used 

TABLE  3 

for  constructing  Histogram 

Category 

Emissivity  Range 

(Bin  Width) 

1 

1.0  -  .875 

2 

.875  -  .75 

3 

.75  -  .25 

4 

.25  -  0 

By  themselves,  histograms  like  Table  3  do  not  give  direct  insight 
into  reflectance  error,  although  knowledge  of  how  many  computed  emissivi- 
ties  fall  into  categories  1  and  2,  for  example,  may  provide  intuitive 
insight  into  the  level  of  error  present  due  to  reflections.  However,  by 
applying  equation  7  in  a  meaningful  manner,  one  may  derive  direct  error 
percentages. 

One  way  of  calculating  meaningful  error  percentages  is  to  assume  a 
worst  case  analysis  and  apply  equation  7  to  each  category  in  the  histo¬ 
gram,  then  multiply  the  percentage  of  emissivities  that  fall  into  each 
category  by  the  error  derived  for  each.  This  results  in  a  worst  case 
weighted  average  error  for  the  view  angle  that  derived  the  histogram.  In 
eguation  form,  this  error  is 


N.*e.  X  100 


(8) 
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where: 

N^.  =  decimal  fraction  of  emissi fives  computed  in  category  i, 

e.  =  error  calculated  using  equation  7  and  lowest  emissivity 
in  category  i  (e.g.,  0.875  for  category  1), 

E  =  worst  case  weighted  average  error  for  the  view  angle  used 
to  derive  the  N^.  in  percent  (%). 

Figure  3  summarizes  the  procedure  used  to  calculate  error  curves. 

3.3  THE  ERROR  CURVES  DERIVED  FROM  THE  PROCEDURE 

Since  looking  at  a  broad  spectrum  of  targets  and  view  angles  would  he 
expensive,  even  using  simulation,  it  was  decided  that  a  single  target,  the 
T72,  and  a  limited  set  of  histograms  based  on  view  angles  around  the  tar¬ 
get  would  be  derived.  These  histograms  were  then  used  as  the  histogram  N^. 
for  equation  8  above.  From  this  limited  data  set,  sufficient  insight  is 
derived  to  draw  further  conclusions. 

Three  azimuths  were  chosen  as  being  representative  of  typical  engage¬ 
ment  scenarios  (which  is  the  ultimate  aim  of  the  total  simulation  task). 
Two  head-on  azimuths  and  an  off  angle  azimuth  were  selected  as  representa¬ 
tive.  Five  elevations  in  all  were  selected.  Three  elevations  at  low 
grazing  angles  were  selected  since  it  was  expected  that  the  greatest  error 
would  be  concentrated  at  these  angles.  High  error  was  expected  at  these 
low  grazing  angles  because,  as  mentioned  earlier,  the  military  paint  simu¬ 
lated  exhibits  specular  behavior  at  low  grazing  angles.  Figure  4  shows 
the  three  azimuths  studied  and  Figure  5  shows  the  five  elevations  studied. 
Tables  4,  5,  and  6  show  the  worst  case  weighted  average  error  computed  for 
azimuth  =  45*,  0*,  and  180*,  respectively,  at  the  elevation  angles  listed 

and  for  the  values  of  reflectance  ratio,  R^,  discussed  in  section  2.0. 

r 

Figure  6  graphs  Tables  4,  5,  and  6. 

Appendix  B  shows  the  actual  histograms  of  emissivities  computed  which 
created  Tables  4  through  6. 
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FIGURE  3.  SUMMARY  OF  PROCEDURE  USED  TO  CALCULATE 
WORST  CASE  ERROR  CURVES. 
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FIGURE  4.  SELECTED  AZIMUTH  ANGLES. 
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Weighted  error  vs.  elevation  angle  for  azimuth  =  45. 


Weighted  error  vs.  elevation  angle  for  azimuth  =  0. 


Weighted  error  vs.  elevation  angle  for  azimuth  =180 


Rr  =  0.0 
Rr  =  0.5 
Rr  =  .667 
Rr  =  1.0 


FIGURE  6. 


ERROR  CURVES  BASED  ON  WEIGHTED  AVERAGE  ERROR  COMPUTED  IN 
TABLES  4,  5,  AND  5. 
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TABLE  4 

Worst  case  weighted  average  Error,  E^, 

for 

Azimuth  =  45° 

El evati on 

Rr=1.0 

Rj,=  .667 

Rr=.5 

Rr=o.o 

0.0 

5.0 

3.18 

5.8 

16.6 

15.0 

5.0 

5.04 

9.1 

23.1 

30.0 

5.0 

2.16 

4.5 

14.0 

60.0 

5.0 

1.39 

2.58 

10.16 

90.0 

5.0 

1.26 

2.19 

9.37 

TABLE  5 

Worst  case 

weighted 

average  Error, 

E„,  for 

Azimuth  =  0.0 

Elevation 

Rr=1.0 

R^=.667 

Rr=-5 

Rr=o.o 

0.0 

5.0 

2.66 

5.5 

13.77 

15.0 

5.0 

6.19 

11.0 

24.7 

30.0 

5.0 

3.75 

6.5 

16.9 

60.0 

5.0 

1.315 

3.4 

9.6 

90.0 

5.0 

1.26 

2.19 

9.37 

TABLE  6 

Worst  case  weighted  average  Error,  E  ,  for  Azimuth  =  180 

W 

Elevation 

Rr=l-0 

R^=.667 

Rr=.5 

R^=0.0 

0.0 

5.0 

1.67 

3.2 

11.4 

15.0 

5.0 

8.59 

14.7 

34.4 

30.0 

5.0 

3.93 

8.11 

21.7 

50.0 

5.0 

1.24 

2.19 

9.33 

90.0 

5.0 

1.26 

2.19 

9.37 
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3.4  INTERPRETATION  OF  THE  ERROR  CALCULATED 

The  curves  in  Figure  6  establish  a  trend  in  error  that  is  view  angle 
dependent.  The  peak  in  error  occurs  at  15*  elevation  for  all  azimuth 
values  and  drops  to  its  lowest  value  at  90*  elevation.  The  reason  for 
this  elevation  angle  dependency  is  due  to  the  specular  behavior  of  the 
emissivity  used. 

Interpretation  needs  to  be  applied  to  the  curves  presented,  since  in 
all  cases,  there  is  error  that  peaks  above  20%  and  in  one  case  peaks  at 
35%.  The  largest  error  is  caused  by  the  curve  computed  for  the  relatively 
uncommon  reflectance  ratio,  R^  =  0.0.  The  other  reflectance  ratios,  0.5, 
0.667,  and  1.0,  are  more  representative  of  realistic  scenarios.  For  these 
other  cases,  the  error  is  mostly  below  10%,  and  peaks  at  15%.  For  eleva¬ 
tion  angles  >  45*,  the  computed  error  is  less  than  5%. 

The  curves  in  Figure  6  represent  error  bounds,  since  no  actual 
vehicle  would  consist  wholly  of  any  one  of  the  reflectance  ratios  chosen. 
In  reality,  the  author  speculates  that  an  average  reflectance  ratio  for  a 
realistic  situation  would  be  somewhere  between  0.667  and  1.0,  based  on  the 
observation  that  many  parts  of  a  vehicle  tend  to  be  near  air  temperature. 

In  summary,  the  following  points  can  be  made: 

•  The  highest  error  occurs  for  a  reflectance  ratio,  R^  =  0.0,  that 
occurs  infrequently  on  a  vehicle  such  as  hot  exhaust  stacks  and 
other  very  hot  engine  parts. 

•  Low  error  was  calculated  for  scenarios  that  might  be  typical  of 
steep  elevation  angles. 

•  Excluding  the  reflectance  case  of  R^  =  0.0,  most  errors  computed 
were  less  than  10%,  which  are  on  the  order  of  6*K  apparent 
temperature  at  300*K. 
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4.0 

SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


Summary 

As  stated  in  section  1.0,  the  purpose  of  this  study  is  to  determine 
whether  it  is  necessary  to  model  reflections  in  the  current  FSTC  method  of 
simulating  IR  signatures.  To  gain  more  insight  the  following  steps  were 
performed: 

1)  A  sensitivity  analysis  approach  was  chosen  for  the  study. 

2)  An  equation  isolating  the  major  components  of  an  IR  signature  was 
derived. 

3)  The  equation  was  used  to  determine  error  present  in  a  simulation 
when  emissivity  was  estimated  to  be  always  0.95.  This  was  accom¬ 
plished  by  using  a  simulation  of  a  complicated  tank  applied  with 
actual,  measured  emissivity  that  had  specular  characteristics. 

4)  A  weighted  average  error  figure  was  derived  and  applied  to  the 
results  of  histogramming  emissivities  for  the  simulated  tank. 

5)  It  was  determined  that  in  some  cases,  error  could  be  as  high  as 
35%  on  an  absolute  scale.  However,  most  of  the  errors  calculated 
in  section  3.0  were  found  to  be  less  than  10%. 

Conclusions 

•  It  was  stated  in  Technical  Memo  #1  that  temperature  is  inferred 
correctly  in  most  of  the  situations  that  FSTC  will  collect  data.  However, 
the  analysis  presented  in  section  2.3  contradicts  this  inference.  There¬ 
fore,  further  study  using  an  approach  similar  to  the  one  presented  must  be 
performed  for  determining  error  in  inferred  temperature.  Unfortunately, 
the  question  of  including  or  not  including  the  effects  of  reflections  when 
inferring  temperatures  from  collected  data  remains  unclear. 

•  It  was  shown  in  section  2.4  that  both  mi s-estimating  emissivity  at 
a  high  value  (~0.95)  and  neglecting  reflections  results  in  compensating 
errors,  if  the  reflectance  ratio,  R^,  is  -1.0.  This  interpretation  is 
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shown  clearly  in  Figure  6  by  the  fact  that  the  greatest  error  is  computed 

for  the  curve  when  =  0. 

r 

Pi scussion 

It  is  inappropriate  at  this  stage  for  ERIM  to  make  value  judgments 
concerning  the  error  in  FSTC's  current  simulation  method,  since  ERIM  is 
unaware  intimately  of  the  needs  and  uses  of  FSTC's  simulation  product. 
However,  it  can  be  stated  that  a  10%  error  is  equivalent  to  about  a  6*K 
apparent  temperature  error  [3].  The  question  of  whether  this  is  important 
or  not  depends  on  the  sensitivity  of  the  IR  systems  being  designed  using 
FSTC's  simulation  products.  Another  key  issue  is  the  fact  that  most  IR 
seeker  systems  are  AC-coupled  devices  that  measure  the  contrast  between  a 
target  and  the  background.  The  analysis  presented  in  this  study  ignores 
this  consideration. 

It  is  possible  to  talk  in  general  terms  about  the  sensitivity  of  IR 
seeker  systems  evaluated  by  ERIM  in  the  past,  as  in  reference  [2]. 

Methods  have  been  used  to  evaluate  IR  system  performance,  usually  defining 
a  probability  of  detection  or  lock-on  of  a  seeker  system  given  range  and 
target-to-background  signal  contrast.  In  some  cases,  small  changes  in 
temperature  contrast  can  result  in  large  differences  in  projected  detec¬ 
tion  ranges. 

Figure  7  shows  curves  relating  lock-on  probability  to  range  at  three 
fixed  temperature  contrasts  for  the  Maverick  seeker  [2].  It  is  striking 
to  the  authors  the  difference  in  these  curves  given  a  single  degree  of 
change  in  temperature  contrast.  If  judgments  are  made  concerning  IR  sys¬ 
tems  with  this  level  of  sensitivity  using  data  that  is  6*  in  error,  the 
judgments  may  be  wrong.  If,  on  the  other  hand,  a  typical  seeker  has  a 
dynamic  range  of  200*K,  6*  of  error  may  not  affect  judgment  at  all. 

Recommendations 

Following  is  a  list  of  recommendations  that  are  based  on  the  conclu¬ 
sions  of  this  report: 
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•  The  sensitivity  of  error  due  to  incorrectly  inferring  temperature 
needs  to  be  studied  in  more  depth.  A  procedure  similar  to  the  one 
used  in  this  report  should  be  used,  as  the  case  study  method  used 
before  is  too  limiting. 

•  A  decision  by  FSTC  personnel  needs  to  be  made  concerning  whether 
the  errors  reported  in  section  3  are  significant  to  their 
customers. 

•  Based  on  the  decision  above,  corrective  action  can  be  planned  for 
by  ERIM  and  FSTC  (i.e.,  including  reflections  and  specular  emis- 
sivity  in  TOTSIG). 
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APPENDIX  A 

DERIVATION  OF  ERROR  EQUATION 

A.l  INTRODUCTION 

In  this  appendix  we  derive  equation  (1)  of  section  2.0,  which  speci¬ 
fies  the  total  error  in  a  simulated  target  radiance.  (Note:  When  we 
speak  of  the  radiance  or  temperature  of  "the  target,"  we  are  actually 
referring  to  some  target  component  or  facet  having  uniform  thermodynamic 
and  optical  properties.)  The  derivation  is  based  on  a  number  of  assump¬ 
tions  regarding  how  the  radiance  simulation  is  performed.  Specifically, 
the  simulation  is  assumed  to  consist  of  the  following  three  steps: 

Step  1.  The  actual  target  radiance  in  some  wavelength  band,  AXj,  is 

measured  by  a  calibrated  radiometer  (e.g.,  a  calibrated  FLIP)  at 
some  view  angle  9^; 

Step  2.  Based  on  the  radiance  measured  in  Step  1,  an  actual  (i.e., 
thermodynamic)  target  temperature  is  calculated; 

Step  3.  The  temperature  from  Step  2  is  then  used  in  a  radiometric  model 
to  predict  target  radiance  at  some  other  view  angle  d^,  and/or 
within  some  other  wavelength  band  AX^. 

A. 2  DEFINITIONS  AND  NOTATIONS 

•  Planck  spectral  radiance  from  a  blackbody,  Lgg(X,  T): 

^  J  C2/XT  J 


C^,  C2  =  constants, 

X  =  wavelength, 

T  =  temperature. 
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In-band  radiance  within  wavelength  band  AX^. ,  L^.  (T) : 
Li(T)  =  f  Lgg(X.  T)  dX 


(A-2) 


•  Differential  in-band  radiance  w.r.t.  temperature, 


It  =  i 
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(Note:  Equations  (A-2)  and  (A-3)  are  typically  evaluated  numerically,  as 
are  many  of  the  equations  given  below.) 

•  In-band  directional  emissivity,  e^.  (5,  T) 


J  e(X,  9,  T)  Lgg(X,  T) 


AX. 

eAe.  T)  =  — 


J  Lbb(X,  t) 


(A-4) 


e(X,  9,  T)  =  spectral,  directional  emissivity  at  temperature  T, 

9  =  abbreviated  notation  for  view  direction  {9,  f). 

•  Bi-directional  reflectance  distribution  function  (BRDF),  p'{9^, 

®i'  = 


4^' 

p'l^r'  K’  ^i’  ^0  "  E(^^. 


(A-5) 
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=  irradiance  incident  from  direction  0^,  (cf.,  Figure  A-1) , 

L(0j,,  (p^)  =  radiance  reflected  into  direction  0^,  (p^  (Figure  A-1). 

Note  that  for  opaque  surfaces, 


£(9.  X)  =  1  -  JJ  9,, 


cos  0^  d0^d^^ 


(A-6) 


Here,  e(0,  X)  is  the  directional  emissivity  in  the  direction  (0., 

•  Reflection  term  into  view  direction  0,  RT(0); 

-  General  surface 


RT 


(9)  -  //  9,.  L(Tj(9,,  cos  9,  d9,d^,  (A-7a) 


Ts(0i,  =  hemispherical  distribution  of  apparent  temperature  (e.g.,  sky 
temperature) . 

-  Perfectly  specular  surface 

RT(0)  =  r(0)  L[T3[0.  -  180*]]  (A-7b) 


r(0)  =  Fresnel  reflectivity, 

Ts(0,  p^  -  180*)  =  apparent  temperature  at  the  specular  geometry: 

di  =  0^:  Pi  =  Pr  -  180*. 

-  Perfectly  diffuse  surface 


4^3(9,,  cos  9,  99, d^. 


(A-7c) 


=  £ 


r 
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r 


28 


ERIM 


FIGURE  A-1 .  BRDF  GEOMETRY  (ADAPTED  FROM  REF.  [4]). 
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p  =  hemispherical  reflectivity, 

Edown  "  tot^l  downwelling  irradiance. 

Note  the  lack  of  any  ^-dependence  for  the  diffuse  case. 

•  Relation  between  actual  and  apparent  temperature 


L.(T^PP(fl))  =  e(e)  +  RT.{S) 


(A-8) 


t3PP(^)  =  target  apparent  temperature  when  viewed  at  angle  9. 

Note  that  in  equation  (A-8)  we  have  not  shown  any  temperature  depend¬ 
ence  for  emissivity.  Although  equation  (A-4)  shows  that  this  is  not 
strictly  valid,  it  is  usually  a  good  approximation  over  temperature  ranges 
of  a  few  tens  of  degrees  Celsius.  This  is  indeed  fortunate  when  tempera¬ 
ture  is  the  unknown  to  be  determined. 


A. 3  DERIVATION  OF  EQUATION  (1) 

Armed  with  the  definitions  of  the  previous  section,  we  may  now 
readily  derive  the  equation  for  error  in  simulated  radiance.  The  deriva¬ 
tion  consists  simply  of  writing  equations  that  describe  the  three  steps  of 
the  simulation  procedure  listed  in  section  A.l. 

Step  1.  Measurement  of  Target  Radiance 

For  a  radiometer  having  a  spectral  passband  AXj  and  view  direction 
9^,  the  measured  radiance  from  the  target  is 

*  RT,(4]  .  (A-9) 


Step  2.  Estimation  of  Actual  Target  Temperature 

FSTC's  present  method  of  temperature  estimation  is  based  on  two 
assumptions:  1)  the  reflected  component  of  the  measured  target  radiance 
is  negligible  (i.e.,  RTj(5j)  =  0);  2)  the  target  emissivity  is  assumed  to 
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be  angul arly-i ndependent  and  have  a  constant  value  of  0.95  in  the  LWIR  and 
0.6  in  the  MWIR.  We  shall  generalize  the  second  assumption  by  writing 

. 


With  these  two  assumptions,  the  estimate  of  actual  temperature,  T, 
may  then  be  written  from  equation  (A-9)  as 


T  = 


(A-10) 


where  L"^(  )  denotes  an  inversion  of  the  band  radiance  versus  temperature 
function  defined  in  equation  (A-2) . 

We  may  express  this  temperature  estimate  as  the  sum  of  the  true  tar¬ 
get  temperature  and  some  temperature  error,  AT,  i.e., 


T  = 


+  AT 


(A-11) 


Note  in  equation  (A-10),  that  T,  and  therefore  AT,  depends  on  the  view 
direction  9y 

Step  3.  Radiance  Simulation 

A* 

The  estimated  target  temperature  T  is  used  to  calculate  target 
radiance  at  some  other  view  angle  and/or  some  other  wavelength  region.  In 
the  general  case,  both  differences  may  apply;  we  therefore  denote  the 

A, 

simulated  radiance  as  12(52)*  again  employ  the  FSTC  radiance  model, 

i.e.,  ignore  any  reflected  component  in  the  new  wavelength  region,  we  may 
wri  te 
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4(^2]  "  ^2  4^^^ 


=  €2  +  at] 


(A-12) 


In  equation  (A-12)  62  is  the  estimate  of  target  emissivity  in  the 
simulated  waveband/view  direction. 

The  correct  value  of  target  radiance  is  given  by  the  following 

4(4]  =  4(4)  4^^^"^^^  ^4(4)  • 

The  error  in  simulated  radiance,  AL2,  is  the  difference  between  (A- 
12)  and  (A-13),  i.e., 

A4  =  4(02)  -  4(4) 


=  €2  +  at]  -  12(^2] 


(A-14) 


act 


If  AT  is  not  too  large,  we  may  approximate  L2(T  +  AT)  by  the  first 
two  terms  of  a  Taylor  series  expansion  around  T  ,  i.e., 

(AT)  .  (A-15) 


P4 

aT 

T^ct^ 

L2[t®‘^^  +  at]  =  L2(t*‘^^]  + 


Inserting  equations  (A-15)  and  (A-13)  into  (A-14)  we  obtain 


9  L 

4L,  =  -  ^2(92))  12^^"'*^  *  '2  ar  ''T  -  ' 
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If  we  consider  the  special  case  of  an  opaque  specular  surface,  the 
reflection  term  may  be  replaced  by  (1  -  £2^^)) 

equation  (A-7b) .  (Note  that  for  an  opaque  surface  the  reflectivity  is  one 
minus  the  emissivity.)  Making  this  replacement  in  the  above  equation,  we 
obtain  our  final  result, 

AL  =  [e  -  e)  +  e  —  AT  -  (1  -  )  LfT^]  .  (A-lo) 

Equation  (A-16)  is  identical  with  equation  (1)  of  section  2.0.  The 
subscripts  and  references  to  6^  have  been  dropped  as  it  is  understood  that 
all  quantities  in  (A-16)  are  defined  for  the  simulated  wavelength  region 
and  view  direction.  The  only  exception  is  the  quantity  AT,  which  is 
independent  of  the  simulation.  Using  equations  (A-9),  (A-10),  and  (A-11), 
we  write  the  following  equation  for  the  error  in  estimated  target 
temperature: 

AT  =  T  -  t""^ 


6^(5) 


(A-17) 


As  shown  above,  AT  depends  only  on  quantities  defined  for  the 
waveband/view  direction  from  which  the  target  radiance  was  actually 
measured. 
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APPENDIX  B 

EMISSIVITY  HISTOGRAMS 

B.l  INTRODUCTION 

In  order  to  compute  the  weighted  average  errors  displayed  in  Tables  4 
through  6,  histograms  of  emissivities,  had  to  be  constructed.  This  was 
done  using  the  procedure  outlined  in  section  3.0.  This  Appendix  lists  the 
actual  histograms  derived  for  the  various  view  angles  considered. 

B.2  COMPUTED  HISTOGRAMS 

Following  Table  3  in  section  3.0,  histograms  of  emissivities  were 
computed.  The  following  tables,  B.l,  B.2,  and  B.3,  show  the  computed 
histograms  for  the  azimuth  and  elevation  angles  listed. 

TABLE  B.l 

HISTOGRAM  FOR  AZIMUTH  =0.0 
Emissivity  Bins 

Elevation  Angle  0-0.25  0.25-0.75  0.75-0.875  0.875-1.0 

0  0.018  0.020  0.460  0.490 

15  0.048  0.188  0.300  0.450 

30  0.020  0.040  0.510  0.410 

60  0.004  0.012  0.090  0.880 

90  0.006  0.008  0.070  0.913 
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TABLE  B.2 

IISTOGRAM  FOR  AZIMUTH  = 

45* 

Emissivity 

Bins 

Elevation  Angle 

0-0.25 

0.25-0.75  0 

.75-0.875  0.875-1.0 

0 

0.073 

0.216  0.689 

15 

0.160 

0.260  0,541 

30 

0.008 

0.015 

0.400  0,570 

60 

0.005 

0.009 

0.134  0.850 

90 

0.006 

0.008 

0.070  0.913 

TABLE  B.3 

HISTOGRAM  FOR  AZIMUTH  =  180* 


Emissivity  Bins 

Elevation  Angle  0-0.25  0.25-0.75  0.75-0.875 


03 
70 
30 
0.004 
0.006 


0.020 

0.306 

0.050 

0.009 

0.008 


0.875-1.0 


.768 
.470 
.270 
.896 
0.913 
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FI3URE  a-1.  HISTOGRAM  FOR  EMISSIVITIES  CALCULATED  AT 
ELEVATION  ANGLE  =  0  AND  ALL  THREE  AZIMUTH 
ANGLES. 
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'IGU’E  2-2.  histogram  -OR  EMISSIVITIES  CALCULATED  AT 

ELEVATION  ANGLE  =  15  AND  ALL  THREE  AZIMUTH 
ANGLES. 
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elevation  =  30 
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FIGURE  B-3.  HISTOGRAM  FOR  EMISSIVITIES  CALCULATED  AT 
ELEVATION  ANGLE  =  30  AND  ALL  THREE  AZIMUTH 
ANGLES. 
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FIGURE  3-4.  histogram  FOR  EMISSIVITIES  CALCULATED  AT 

ELEVATION  ANGLE  =  60  AND  ALL  THREE  AilMUTH 
ANGLES. 
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FIGURE  8-5,  HISTOGRAM  FOR  EMISSIVITIES  CALCULATED  AT 
ELEVATION  ANGLE  »  90  AND  ALL  THREE  AZIMUTH 
ANGLES. 


38 


